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PARTITIONING AND PACKING MATHEMATICAL SIMULATION MODELS 


FOR CALCULATION ON PARALLEL COMPUTERS 


Dale J. Arpasl and Edward J. Milner 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

This paper deals with the development of multiprocessor simulations from 
a serial set of ordinary differential equations describing a physical system. 
Degrees of parallelism (i.e., coupling between the equations) and their Impact 
on parallel processing are discussed. The problem of identifying computational 
parallelism within sets of closely coupled equations that require the exchange 
of current values of variables is described. 

A technique is presented for identifying this parallelism and for parti- 
tioning the equations for parallel solution on a multiprocessor. An algorithm 
which packs the equations into a minimum number of processors Is also 
described. The results of the packing algorithm when applied to a turbojet 
engine model are presented In terms of processor utilization. 


INTRODUCTION 

Multiple processors, operating together to solve a single problem, can. 

In many cases, decrease the time of calculation. This Is Important In time- 
critical applications, such as real-time simulation, where this technique can 
provide computational rates unachievable on a single processor or allow the use 
of lower cost hardware to provide the necessary computational capabilities. 

For certain classes of problems It Is possible to configure a network of micro- 
computers to achieve the same throughput rate as a large mainframe computer at 
a lower Initial and ongoing maintenance cost. 

The parallel processing concept has opened new areas of research and 
development In hardware, software and theory. Some efforts sponsored by NASA 
Lewis are described In references 1 to 4. Techniques for developing mathe- 
matical models that can be solved efficiently on parallel processors Is a key 
area of study. The first step In developing these multiprocessor models. Is to 
Identify parallelism within the mathematical formulation of the problem. This 
requires a data flow analysis of the problem's equations and will Identify the 
"critical path" and the minimum achievable calculation time. The next step Is 
to arrange, or "pack" the noncrltical path computations on the minimum number 
of processors so as to make maximum use of the available computing resources. 

This paper presents a method of partitioning equations for multiprocessor 
solution. The method, when applied to models containing Inherent parallelism, 
will produce segmented sets of equations that can be solved In parallel. A 
brief discussion of computational parallelism Is presented In terms of the 
degree of coupling between subsets of the model equations. The Identification 
of parallel paths In closely coupled sets of equations Is then discussed, 
followed by the description of a packing algorithm for assignment of the paths 







to a minimum number of processors. Finally, the results of the algorithm 
applied to the mathematical model of a helicopter engine are used to demon- 
strate the algorithm's capabilities. 


COMPUTATIONAL PARALLELISM 

A mathematical model of a physical system consists of a set of equations 
which describe, to some degree of accuracy, the response of that system to 
external Influences (driving functions) over a limited range of operation. 

This range Is defined In terms of the maximum and minimum values of the driving 
functions and. If time dependent, the maximum frequency or maximum rate of 
change of these functions. Generally, the object of this modeling effort Is 
to provide a simulation of the physical system. 

Whenever a simulation must Interface to the real world (e.g., with control 
hardware), the mathematical model must be computed at a rate sufficient to make 
the simulation look like the real physical process to the real-world compo- 
nents. This Is termed real-time computation. To test a piece of control hard- 
ware, the simulation response to control Inputs must dynamically match the 
response of the physical system. This real-time requirement establishes a max- 
imum allowable calculation time for a digital simulation. The maximum allow- 
able calculation time may be further reduced by Internal numerical stability 
requirements. The maximum calculation time necessary to meet all of these 
real-time requirements Is termed the required update time (AT). The simulation 
will meet all real-time requirements If the mathematical model Is computed 
every AT sec. 

Once the simulation update time has been established, the model must be 
programmed and executed on a computer. Both the efficiency of the programming 
language and the capabilities of the computer hardware will determine the 
achievable calculation time of the model. If the model can be divided Into 
parts which can be computed In parallel, then the efficiency requirements of 
the language and the computational capabilities of the hardware may be relaxed 
by assigning each path to a separate computer and computing these paths In 
parallel (assuming appropriate data transfer takes place between the com- 
puters.) If a model can be computed serially In the required update time on a 
single computer, parallel processing may allow the use of a number of lower 
cost computers to do the job. This could be a more cost effective approach to 
simulation. 

The prerequisite to developing parallel processor simulations Is to be 
able to Identify the parallel computational paths contained In the model. In 
general, a dynamic model can be programmed on a digital computer as a set of 
N computationally sequential equations of the form 

X K Oh) = f kE X m ( 1 h ) , X m ((1 - l)h) u(1h) ] 

where Xj((1h) Is the result of the K*h equation at time 1h. Here, h 
denotes the simulation time step or update Interval of the model calculations. 
The arguments X m (1h) are the current values of the results of preceding 
equations In the model (l.e., m = 1 to K - 1), and X m ((1 - l)h) are the 
past values of the results of all equations In the model. The argument u 
represents values obtained from sources external to the model which are always 
available at the start of the model computation sequence. The functional 
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relationship between X« and Its arguments is represented by f|<. Assuming 
an equation Is an Indivisible computational unit, then the parallelism In the 
model Is determined by the arguments of each equation. That Is, two equations, 
or sets of equations, can be computed In parallel within an update Interval If 
their arguments are Independent of the results of the others computed In that 
Interval. 

For example, a model of the form 

Xi(lh) = f 1 [X 3 ((1 - l)h)] 

X 2 (1h) = f 2 [X 1 ((1h))] 

X 3 (1h) = f 3 [X 2 (( 1h) ) , X 3 ((1 - 1 )h) ] 

contains no parallelism since X 3 ( 1 h) requires X 2 ( 1 h) and X 2 ( 1 h ) requires 
X-|(1h). These calculations must be done serially. However, the model 

Xidh) = f t [X 3 ( ( 1 - l)h)] 

X 2 (1h) = f 2 [u(1h)] 

X 3 ( 1 h ) = f 3 [X] ( 1h) , X 2 (1h), X 3 ((1 - l)h] 
does contain parallelism since X] can be computed at the same time as X 2 . 

Many times, parallelism Is obvious from the physical nature of the system 
which Is modeled. For example, to test a control system designed to balance 
the load on two generators, the generators could be formulated In a model con- 
taining two sets of decoupled equations. That Is, there are no result values 
transferred between the sets. Both, however, provide results to an external 
piece of hardware (the controller). It Is obvious that the simulation of the 
generators can at least be Implemented on two computers operating In parallel 
(fig. 1(a)). Each equation set (generator) may, of course, contain additional 
computational parallelism. 

Another example of physically detectable parallelism occurs when two or 
more sets of equations have significantly different dynamic characteristics. 

For example, suppose a simulation of a single engine aircraft Is desired. The 
mathematical model would consist of a set of engine equations and a set of 
airframe equations. The airframe dynamics are generally an order of magnitude 
slower than the engine dynamics. Therefore, these equation sets can be com- 
puted using different update times. If the update times differ by a factor of 
5 or more (rule of thumb), then results transferred between the airframe and 
engine can be past values without Introducing dynamic errors. In this case, 
the equation sets are considered loosely coupled and each can be placed on a 
separate computer. As shown In figure 1(b), because the airframe equations 
(A) use only results of engine equations (E) computed In a previous calculation 
(update Interval) and vice-versa, these arguments can be computed In parallel. 
As Is the case of decoupled equation sets, loosely coupled equation sets may 
each contain additional computational parallelism. 
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Parallelism due to decoupled or loosely coupled equation sets Is easily 
Identified from the physical nature of the model. A more difficult task Is the 
Identification of parallelism In a set of closely coupled equations, where the 
process dynamics dictate the use of current arguments In solving for equation 
results. For Instance, suppose a model contains the following set of 
equations: 


Xidh) = f 1 [X 5 ((1 - l)h), u(1h)] 

X 2 dh) = f 2 [X 5 (d - l)h)] 

X 3 (1h) = f 3 [X 2 (1h)] 

X 4 (1h) = f 4 [ X-j ( 1 h ) , X 2 (1h) ] 

X 5 ( 1 h) = f 5 [X 3 (1h), X 4 < 1 h) , X 5 ((1 - 1 ) h) ] 

The variable X-| can be computed at the start of the calculation Interval, 
since It Is a function of the past value of X5 and the external variable 
u. X 2 may also be computed at the start of the Interval. However, the 
calculation of X 3 must be delayed until X 2 has been determined, the cal- 
culation of X 4 must be delayed until both X-| and X 2 are determined, and 
the calculation of X5 must be delayed until both X3 and X 4 have been 
determined. As shown In figure 2(a), three computational paths can be Identi- 
fied which can be assigned to three different computers In the simulation. 

Note that "wait states" have been Inserted to Insure the currency of the argu- 
ments. That Is, equation calculation Is delayed until current argument values 
become available. The X3 calculation Is shown to be delayed slightly for 
the transfer of X 2 . The shaded areas (time slots) Indicate the time available 
for result transfer to computer number 1. The calculation of Xi and X3 
can take place anywhere In the time slot. 

The detection of this type of computational parallelism can become burden- 
some when the equation set becomes large. The technique, however, can be auto- 
mated. Related to this problem of partitioning Is the problem of allocation 
(l.e., packing these paths Into a minimum number of computers without extending 
the update time). Figure 2(b) demonstrates packing of the paths defined In 
figure 2(a). Arbitrary calculation times of TX-j , TX 2 , TX3, TX 4 and TX5 have 
been assigned to the equations producing results Xi through X5, respec- 
tively. The time TX-| Includes the time required to obtain the value of u. 
Note In figure 2(b), that, because of the calculation times, the X 2 - X 4 - X5 
path Is critical In that It contains no Idle states. This path, therefore, 
dictates the minimum possible update time (AT = TX2 + TX 4 + TX 5 ). The paths 
X-j and X3 are assigned to separate computers. Packing In. this example Is a 
trivial task, since the X 3 calculation can be moved onto the computer 
number 2 to be calculated during the Idle period. 

In many cases, efficient packing requires shifting equations In their 
time slots. This causes a ripple effect on the time slots of other equations 
which can complicate the packing problem. Because of the nature of the packing 
problem, a unique solution to the development of a packing algorithm does not 
exist. There are many ways to pack most parallel models. The optimum approach 
may be model-dependent. 
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In the following sections partitioning and packing algorithms that have 
been developed at NASA Lewis are discussed. These algorithms were tested with 
a model of a jet engine and the results are presented. 


PARTITIONING 

To begin the discussion of the partitioning algorithm, certain terms 
should be defined. A mathematical model Is a set of equations, written to 
define the characteristics of a physical system to some desired degree of 
accuracy. A program Is a sequential set of digital equations and supporting 
Information (e.g., variable and constant definition) which define the mathe- 
matical model within the constructs of a programming language. A path Is a 
subset of these equations which, because of Interrelationships between argu- 
ments and results, contains no parallelism. Partitioning Is the transformation 
of the program equations Into a number of paths which may be calculated In 
parallel. Packing Is the combination of paths Into a minimum number of pro- 
cessors (computers) which, provide computation of the model within a prescribed 
update Interval. The critical path Is the longest path and the prescribed 
update Interval must be greater than or equal to the calculation time of the 
critical path. 

In this discussion of partitioning It Is assumed that a program Is given. 
That Is, these equations, when executed serially, provide the required results. 
No assumptions are made concerning the parallelism of computational units 
(operations) contained In the program equations. The equation 

x = a*y + b*z 

contains parallelism (l.e., a*y can be calculated In parallel with b*z) which 
will be Ignored since we are concerned with partitioning at the equation level 
and not parsing. For purposes of this discussion, the above equation will be 
considered as 


x = f (a , y, b, z) 

where f Is some single operation. Therefore, equations will be assigned to 
paths In their entirely and not broken up Into more primitive result-argument 
relationships . 

As Indicated In the last section, partitioning requires the establishment 
of result-argument relationships for the serial set of equations In order to 
develop computational paths. It Is also necessary to know the calculation time 
of each equation. The program must be processed to provide this Information. 
For this effort, the result-argument relationships and the calculation time 
Information are outputs of the multiprocessor programming utility RTMPL 
(refs. 2 and 3). The primary function of this utility Is to translate a struc- 
tured program of the mathematical model Into assembly language for the simula- 
tion processor(s) . As an option, the utility also provides Information on the 
result and arguments of each equation, the processor operations necessary to 
obtain the result, and the processor calculation time for each operation. For 
example, the utility-generated Information for the equation 
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X = y + z 


might be as shown In table I. Each equation has a label (programmer or RTMPL 
defined). In this case, S$12 Is used to Indicate that this Information con- 
cerns the twelfth equation In the program. The processor operations to compute 
the equation are: load register R1 with z (requires 8 time units), add 
variable y to R1 (16 time units) and store R1 as the value of variable 
X (8 time units). This type of Information Is generated for each equation In 
the program. 

The first step In the partitioning process Is to convert the utility gen- 
erated Information Into the form needed for partitioning. This form Is shown 
In table II for the close-coupled example In the previous section. To simplify 
the discussion, the equation label has been replaced by the name of the result 
variable. Dependent arguments are those which are the results of previous 
program equations calculations In the update Interval (e.g., X-] Is a depend- 
ent argument of X 4 ). These are the drivers for partitioning since their 
current values are required before the computation sequence can continue. The 
Independent arguments u and X 5 do not affect partitioning since only past 
values are used. The calculation time for each equation Is determined by add- 
ing the calculation times of the given processor operations. For example, the 
calculation time of equation X-| Is determined to be 32 from table,!. 

The time at which an equation can start Is determined by the arguments and 
calculation time of each equation. The first equation of a set only has Inde- 
pendent- arguments and thus, can always start at time 0 (measured from the 
beginning of the calculation update Interval). It can never require results 
from calculations In the current update Interval since none are yet available. 
An equation can end at the time obtained by adding Its calculation time to the 
time It can start. The general formula for obtaining this time Is 

CANSTART( RESULT) = MAX(CANEND(ARG 1), CANEND(ARG 2), ..., 0) 

CANEND( RESULT) = CANSTART( RESULT) + CALCTIHE( RESULT) 

where ARG1 Is the first dependent argument, etc. This formula Is applied 
sequentially to each equation In the program (see table II). 

Once these attributes have been established for each program equation, the 
Identification of computational paths contained In the program can begin. The 
algorithm used for path Identification Is shown In figure 3. Its purpose Is to 
Identify all sequences of equations which contain no parallelism and which must 
be computed serially. These paths are organized Into a linked list called 
PATHLIST. The paths In PATHLIST are ordered In terms of decreasing path cal- 
culation time. Therefore, the first path In PATHLIST Is the critical path. 

To form a path, the algorithm selects the program equation, having the maximum 
CANEND time, and which has not already been assigned to a path. This Is the 
result equation of the path. The next equation selected Is the one which 
produces a result used as a dependent argument of the result equation. If more 
than one equation result Is used as a dependent argument, then the one with 
maximum CANEND time Is selected. The selected equation then Is Inserted In 
front of the result equation In the path. The path formation process continues 
until an equation Is Inserted which has no dependent arguments equations which 
are not already assigned to a path. Paths are formed until all program 
equations have been assigned. 
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Partitioning has been discussed In terms of equations that produce values 
of variables. Often, mathematical models contain statements that do not pro- 
duce values. Two common examples are conditional statements (e.g., IF ... 

THEN ... ELSE) and command statements (e.g., I/O operations). The calculation 
time of such a statement must be combined with a preceedlng or following equa- 
tion. This could Impose limitations on program structure and Is a subject for 
future study. 


PACKING 

The partitioning process produces a number of paths consisting of equa- 
tions which must be computed serially and a table of Information on each equa- 
tion (described In table II). The final task In the process of formulating a 
multiprocessor model Is to pack the paths for assignment to a minimum number of 
processors. The first path In our list has the largest calculation time due to 
the partitioning algorithm. This Is called the critical path and Its calcula- 
tion time Is the minimum time In which the model can be computed no matter how 
many processors are used. The number of paths Identified through partitioning 
Is usually greater than the minimum number of processors necessary. 

The minimum number of processers necessary to Implement a multiprocessor 
simulation depends on how fast the simulation must be computed. This update 
time must be specified prior to packing. The simulation time step h Is 
usually based on stability and dynamic accuracy requirements. For real-time 
applications, the update time AT must be equal to h. The update time also 
specifies when the computations must end. The first step In packing Is to 
determine when each equation must end using the specified update time. By 
doing so the last column of Information Is added to table II. 

To determine when an equation must end, we begin with the state variables 
(defined here as those variables whose current values are not used as arguments 
In the model, but appear as results of model equations). The state variable 
computations will be the last computations performed, and thus must will end 
at the prescribed update Intervals. The calculation of equations which are 
dependent arguments of these variables must end no later than the time at which 
the state variable calculations must start. The times when subsequent equa- 
tions In the result/argument string must end Is similarly determined. Since a 
variable can be used as a dependent argument In more than one equation, care 
must be taken that the earliest time, arrived at after all paths are analyzed, 
Is used to specify when that equation must end. The numbers given In table II 
were arrived at by specifying an update time equal to the time at which the 
state variable equation X5 can end. Note that both X3 and X4 use X2 as 
a dependent argument. The time at which X2 must end, as determined from X3 
requirements. Is 

MUSTEND(X 2 ) = 112 - 32 = 80 

and, from X4, Is 


HUSTEND(X 2 ) = 112 - 48 = 64 


The minimum Is selected. 
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We now have determined when an equation can start, can end, and must end. 
These are termed equation attributes. Since the paths are serial they can also 
be assigned these attributes: A path can start when Its first equation can 

start, a path can end when Its last equation can end, a path must end when Its 
last equation must end, and additionally, the calculation time of a path Is the 
summation of the calculation times of Its equations. This Is sufficient 
Information to pack the paths. 

The solution to the packing problem Is not unique In that many arrange- 
ments of paths In processors can result In a satisfactory solution. The 
requirements placed on a general packing algorithm are not strict, however, 
from the point of view of efficiency of processor utilization. For example. It 
does not matter If the last processor that Is packed performs Its calculations 
In 5 or 95 percent of the update time If there Is Insufficient time to calculate 
all of Its equations on the other processors used In the simulation. From a 
processor utilization point of view, both packing arrangements are satis- 
factory since the same number of processors are used In the simulation. 

The packing algorithm, shown In figure 4, was designed to achieve the 
minimum number of processors. Other requirements which may be Imposed, such 
as memory size limitations and Inter-processor data transfer limitations were 
not Imposed on the algorithm. 

As Input, the algorithm requires: 1) that all paths be specified In a 

linked list called PATHLIST In order of decreasing calculation time; 2) that 
the required update time of the simulation, AT, Is specified and that the 
attributes of each equation and path (CANSTART, CANEND, MUSTEND, CALCTIME) 
have been determined as described above. 

The packing algorithm creates processors as needed and Inserts paths from 
PATHLIST according to a hierarchy of relationships between existing equations 
In a processor and the equations In the unpacked paths. When a processor Is 
created, the path with the longest calculation time In PATHLIST Is Inserted. 
Next, the paths which are related to paths already In a processor are tested 
to see If they fit (see discussion of TESTFIT algorithm below). If so, they 
are Inserted, If not, they are placed In a carry-over list. 

Then, paths In PATHLIST which are unrelated to the equations In the pro- 
cessor are tested. If one of these Is Inserted, unrelatlonal testing Is ended 
and relational testing begins again. When no other paths can be Inserted Into 
a processor, another processor Is created. This process continues until all 
paths In PATHLIST are Inserted Into a processor. 

Relational testing Is prioritized. All unpacked paths which provide 
critical arguments are tested first. (A path Is considered to provide a 
critical argument If the result of the last equation In the path (EL) Is an 
argument of a processor equation (EP) and 

MUSTEND (EL) = MUSTEND (EP) - CALCTIME (EP) 

Next, other related paths are tested. Then paths In the carry over list (which 
was formed from paths which were related to equations packed Into previously 
formed processors, but not yet packed) are tested. 
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Paths are tested for Insertion on an equation by equation basis using the 
test fit algorithm shown In figure 5. First the attributes (CANSTART, CANEND, 
CALCTIME, and MUSTEND) of all program equations are saved. This Is necessary 
because Inserting an equation Into a processor can cause a ripple effect on the 
attributes of other equations. If the whole path does not fit, any equation of 
the path, Inserted Into the processor, must be removed and the attributes of 
affected equations restored. 

The ripple effect Is Illustrated In figure 6. Assume a processor contains 
two equations (A and B) and that It has been determined that equation (C) can 
be Inserted between them. The calculation time of each equation Is shown as 
the shaded areas. For packing purposes, the calculation of each equation can 
take place anytime between Its CANSTART time and Its MUSTEND time. 

The space available for C Is the difference between the time at which B 
must end and A can end minus the calculation time of B. . Equation C will be 
Inserted to start directly after A can end. The calculation of B will be 
delayed until C can end. Note that the difference between the MUSTEND and 
CANEND times of A and B have been reduced to zero by the positioning of C, and 
that the time difference for C has been reduced. The primary Impact of these 
changes Is to reduce the space In the processor available for packing other 
paths. There Is also a secondary Impact of equal Importance. By Increasing 
the times at which C and B can end, any unpacked equations which use these 
equations as arguments have their starting times delayed. Similarly, by reduc- 
ing the MUSTEND times of A and C, the MUSTEND times of any unpacked argu- 
ments of these equations are moved up. These effects tend to reduce the slot 
sizes of unpacked equations restricting the range of time Into which they can 
be packed Into a processor. Also, these ripple effects may Introduce computa- 
tional gaps within unpacked paths. 

After the attributes are saved (fig. 5), the path equations are ordered 
In terms of decreasing CANEND times for Insertion testing. That Is, the 
latest equation will be tested first and the earliest last. 

The processor equations are arranged In sequential order where EP(1) Is 
the earliest equation and EP(n) Is the latest equation. Testing to determine 
If a path equation ( E( 1 ) ) can be Inserted Into the processor Involves the 
Identification of all slots between any two processor equations (EP(j - 1), 
EP(j)) where the equation fits. The processor end points ( 1 . e . , EP(j) = EP(1) 
and EP(j - 1) = EP(n)) must also be considered. Because of the argument and 
result relationships between E(1) and the processor equations It Is required 
that the range of processor equations be limited for testing purposes. Let the 
end points of the range be designated by EPE and EPL (the earliest and 
latest processor equations respectively, before which E(1) may be Inserted). 
This range Is established as follows: If E(1) Is an argument of a processor 

equation, then EPL Is the earliest processor equation of which E(1) Is an 
argument (EPE = EP(1)); If any processor equation Is an argument of E(1), 
then EPE Is the one following the latest of these and EPL Is the last pro- 
cessor equation plus one (end point); If E(1) Is unrelated to any processor 
equation, then EPE = EP(1) and EPL Is the last processor equation plus 
one. 


Once the range of testing has been established, all slots within that 
range are tested to determine If E(1) fits. The fit criterion Is as follows: 
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1. If EP( j) = EP(1 ) then CANEND* (E(1)) = CALCTIME (E(1)) else CANEND* 
( E ( 1 ) ) = CANEND (EP( j - 1)) + CALCTIME (E(1) ) ; 

2. If CANEND* ( E( 1 ) ) < CANEND ( E( i ) ) then CANEND* ( E( 1 ) ) = CANEND ( E( 1 ) ) ; 

3. If EP( j - 1) = EP(n) then MUSTEND* ( E( 1 ) ) = AT else 
MUSTEND* ( E( 1 ) ) = MUSTEND (EP(j)) - CALCTIME (EP( j ) ) ; 

4. If MUSTEND* ( E( 1 ) ) > MUSTEND ( E( 1 ) ) then 
MUSTEND* ( E( 1 ) ) = MUSTEND (E(1)) 

5. IF [[CANEND* ( E( 1 ) ) < = MUSTEND (E(1))] and 
[CANEND* ( E( 1 ) ) < = MUSTEND* ( E( 1 ) ) ] ] then E(1) fits. 

The asterisk Indicates the attributes of E(1) If It were Inserted Into the 
slot between EP(j - 1) and EP(j). 

If It Is established that E(1) can fit In more than one slot, the test- 
fit algorithm proceeds to select the slot Into which E(1) best fits. The 
best fit criterion Is as follows: 

If a slot exists such that 

CANEND* ( E ( 1 ) ) - CALCTIME (E( 1 ) ) - CANEND (EP(j - 1)) = 0 
then this slot Is selected. Otherwise, the latest slot which maximizes 

[MUSTEND* ( E( 1 ) ) - CANEND* (E(1 ) ) ] 

Is selected. 

This criterion provides for efficient packing by eliminating processor 
Idle time If possible, and If not, then the ripple effect from Insertion Is 
minimized. 

Once a slot Is selected, the equation Is Inserted and the attributes of 
all program equations are updated to reflect the Insertion. If any path equa- 
tion cannot be Inserted Into the processor, path equations which have already 
been Inserted are removed from the processor, the original attributes are 
restored to the program equations and the Test Fit algorithm ends. 


RESULTS 

The packing algorithm was programmed In Pascal, along with the partition- 
ing algorithm. It was then tested on a helicopter engine model. The appendix 
contains a complete listing and description of the output. The results, In 
terms of percent processor utilization, are presented In table III. This first 
column Is the update time AT specified prior to packing. It Is given In 
terms of machine cycles. The second column gives the number of processors 
required for packing. The remaining columns show the percent utilization of 
the update time In calculating each processor's assigned equations. The first 
specified update time (5666 cycles) was the minimum possible time and corre- 
sponds to the critical path calculation time. In this case, four processors 
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were required. The second update time selection (10 000 cycles) required two 
processors. Note that the percent utilization of the last processor In each 
case exceeds the summation of time available on the other processor(s) . The 
algorithm, therefore, functioned satisfactorily. 

Since the packing algorithm does not account for data transfer times, It 
Is possible that the available time between when a variable Is computed on one 
processor and when Its value Is required for computation on another processor 
will be less than the time required to transfer the variable between the pro- 
cessors. The effect of this will be to Increase the effective calculation 
time of the packed simulation, and therefore, to Increase the minimum achiev- 
able update time. Data transfer effects may be significant for multiprocessor 
systems with Inefficient data transfer mechanisms or for simulations that 
require large volumes of data transfer between processors. Future work In the 
development of a packing algorithm should Include a study of the effects of 
data transfer. While these effects will Increase the critical path time, and 
therefore, the minimum update time, proper consideration of data transfer will 
minimize this Increase and provide for more efficient packing. 


CONCLUDING REMARKS 

The algorithms and considerations presented for partitioning and packing 
mathematical models for calculation on parallel processors has simplified the 
development of multiprocessor simulations at NASA Lewis. Evaluation of the 
packing algorithm will continue as other multiprocessor simulations are devel- 
oped. Work on a completely automated programming package Is In progress, which 
will take a structured serial statement of a mathematical model, detect Its 
parallelism, and provide load modules for the required number of processors. 

The authors welcome discussions of the techniques presented In the paper, 
related techniques, and developments In the many other aspects of multiproces- 
sor simulation. 
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APPENDIX 


The following listing Is a result of partitioning and packing the program 
equations describing a helicopter engine. An update Interval equal to the 
critical path time was specified. All times are given In terms of machine 
cycle time of the Motorola 68000 microprocessor. 

Listing pages 1 to 3 list the program equations In terms of result vari- 
ables. For each equation, the calculation time (EXTIME) and the CANSTART 
and MUSTEND attributes are given. The "RELATED EQUATION" column Identifies 
the relationships between the equation and other equations. If the equation Is 
an argument of other equations It Is so Indicated by an "UB" entry. Constant 
arguments of the equation are designated by "CN" followed by the name of the 
constant. An "SV" entry designates those equations whose past values are argu- 
ments of the equation. Current value or dependent argument equations are 
Indicated by a "RQ" entry. 

The paths generated by the partitioning algorithm are shown at the top of 
listing page 4. Nineteen paths were found. The attributes of each path are 
given, followed by the equations contained In each path. Path number 1 Is the 
critical path. 

The packing sequence beginning In the middle of listing page 4 and ending 
on page 11 Is an optional diagnostic listing which can be selected by the user. 
It details the operation of the packing algorithm In processing the paths. 

Each step Is preceded by asterisks. The step Is then defined and the operable 
equations Identified. The number preceding the equation list Is the appro- 
priate path number. If applicable. The number zero Implies that the equations 
are contained In a temporary working list. The last comment on listing page 11 
Indicates whether or not the processors were successfully packed. If the 
smallest path (processor) execution time Is less than the unused time on the 
other paths (processors) then the packing Is successful. 

Listing pages 12 to 14 define the equations packed Into each processor. 

The processors are Identified by number (In this case 1 to 4) and the percent 
utilization of the prescribed update Interval Is given. The calculation time 
(CALC), the time (relative to the start of the update Interval) at which the 
equation starts and ends Its computation, and the time computation must end 
(MAX) In order to meet data transfer requirements Is listed. These require- 
ments are also listed for each equation. The processors ("PR") and their equa- 
tions which use an equation as a dependent argument are shown In the "SENT TO" 
column. The time at which the value Is required to arrive at the destination 
processor (REQ) is also provided. The “NEEDS" column Indicates those proces- 
sors and equations which supply dependent argument values to an equation. The 
time at which these values are available (AVAIL) Is also Included. 
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PRTRTMPL LISTING 


M 6800 Q 


02/19/85 


1^3 f U16 


PAGE 1 


EQUATION 

EXTIME 

CANSTRT 

MUSTEND 

RE!... 

ATED EQUATIONS 

AEI...2 

218 

0 

3856 

UB ; 

WA2 





CNJ 

P2 

RTTH2 

770 

0 

77 0 

!JB J 

WA2 





cnj 

T2 

WF 

2 IT 

0 

2316 

UBJ 

FART1 





CN J 

WFPH 

P3 

2*1 2 

0 

1076 

UBJ 

PS3 





sin 

T3--WS8 

PS3 

122 

2T2 

1198 

Rt:n 

P3 

PS3Q2 

20 T 

36T 

IT-02 

Ren 

PS3 





UBJ 

WA2C-T25Q2 





CN J 

P2 

T3Q2 

598 

563 

T832 

RQJ 

PS3Q2 





GNJ 

XPRC1-NFRC1-ZTRC 

T3C 

1ST 

1 166 

T966 

RQJ 

T3Q2 





CN * 

T2 

NGC 

T3T 

770 

12 OT 

RQJ 

RTTH2 





sin 

NG 

PCNGC 

198 

120*1 

IT 02 

RQ J 

NGC 





UBJ 

B1 

WA2C 

2136 

1*102 

3538 

RQJ 

PS3Q2-PCNGC 





UBJ 

WXG2-WA2 





CN J 

XPRC—YFNGC— NPRC—ZW2C 

WA2 

292 

3533 

T1T8 

RQJ 

MA2C-DEL2-RTTH2 





UBJ 

W B 2 5 - W A 3 -- TO R Q C— P T 5 D T~ W 1! 

PI 

598 

IT 02 

T1T8 

RQJ 

PCNGC 





UBJ 

WB25 





CN J 

XPNGC— NPNGC— ZB1 

B2 

610 

3538 

TITS 

RQJ 

WA2C 





CN J 

XW2CB-NW2CB-ZB2 

WB25 

I 92 

T1T3 

•T270 

RQJ 

B1-B2-WA2 

NAS 

30 

T270 

T30 0 

RQJ 

MA2-MB25 





UBJ 

WS3DT 

HXQ2 

598 

3533 

T990 

RQJ 

WA2C 





UBJ 

MB3 





CN J 

XW2C-NW2C-ZWXQ2 

m3 

130 

T 1 36 

5330 

RQJ 

WXQ2-WA2 





UBJ 

WS3DT 

WAS 1 

970 

2T2 

2316 

RQJ 

P3 





UBJ 

PT1DT-W83DT 





SV J 

PT1-WS3 

H3 

IT 2 

0 

25T0 

UBJ 

TORQC-HT1 





SV J 

T3 

FART1 

22T 

1212 

25T0 

RQJ 

WF-WA31 

HT1 

336 

IT 36 

2876 

RQJ 

FART1--H3 

TT 1 

122 

1772 


RQJ 

HT1 





UB J 

PT1DT 

TI-ITAT1 

122 

1 89T 

3120 

RQJ 

TT1 





UBJ 

WTl 

NT 1 

838 

2016 

T50T 

RQJ 

THTAT1 

. 




UBJ 

PT5DT -TGRQT 1 





sin 

PT1 

'"•RTSQ:!. 

230 

0 

'•> 

/•. 

SV J 

PT1--PT5 

HQTHT 

598 

230 

3120 

r\\A * 

PRT5Q1 
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M 680 00 


02/19/as 14:34: 16 


PRTRTMPL 

LISTING 

.-141 

118 

T0RQ41 

390 

T25Q2 

134 

T25 

lie 

1-125 

114 

H2 

126 

TORQC 

594 

H44 

54 

1-145 

114 

T45 

122 

PR49Q5 

230 

W45C 

598 

T HTA45 

122 

W45 

854 

D1-IQTI-I5 

610 

DI-145 

118 

T0RQ45 

390 

1-149 

X. JO 

NGDT 

566 

NPDT 

566 

P41DT 

312 

P45DT 

470 

T3DT 

494 

WS3DT 

130 

NG 

206 

NP 

206 

•LI. 

206 




UB J 

cn: 

2016 

3238 

RQJ 

ub: 

2854 

4894 

RQ % 

ue: 

sv: 

568 

4 0 68 

RQJ 

702 

4186 

RQJ 

cn: 

820 

4300 

RQJ 

ub: 

0 

4300 

ub: 

cn: 

4300 - 

4894 

rq: 
s v: 

2134 

3292 

rq: 

2188 

3406 

RQJ 

ub: 

2302 

3528 

rq: 

ub: 

0 

3052 

ub: 

cn: 

sv: 

230 

3650 

RQJ 

ub: 

cn: 

2424 

3650 

rq: 

ub: 

2546 

4504 

RQJ 

ub: 
sv : 

230 

4386 

RQJ 

ub: 

cn: 

2546 

4504 

rq: 

ub: 

340 0 

4894 

rq: 

sv: 

2664 

5666 

RQJ 

4894 

5460 

RQJ 

3790 

5460 

rq: 

cn: 

2854 

5460 

rq: 

4136 

5460 

RQJ 

130 0 

5460 

rq : 
su: 

4300 

5460 

RQJ 

5460 

5666 

rq : 
cn: 
svj 

4356 

5666 

rq : 
cn: 
s 1 , 1 : 

3166 

566 A 

RQJ 

cn: 


PAGE 2 

Dll '11 

XPRGT-NFRGT-ZDHGTQ 

DI-luTl-l 4-THTA4 1 

TORQ41 

W41-DH41 

NGDT 

NG 

PS3Q2 

T25Q2 

T75 

TORQC 

torch: 

T2 

H25-WB25-H2-WA2-H3-WA3 

NG 

H41-DH4! 

1-144 

H49 

1-145 

P45DT 

W45C 

P49 

P45 

PR49Q5 

W45 

XPRPT-NPRPT-ZW45C 
T45 
D 1-145 

T HT A 4 5 - W 4 5C 

P45DT 

P45 

FR49G5 

DH45 

XPRPT-NPRPT— ZDHPTQ 

DHQTH5--THTA45 

T0RQ45 

W45-DI-I45 

NP 

H45-DII45 

T0RQ41-TQRQC 

TDRR45 

TORQLD 

W41-WA31-T41 

WXQ2-WA2-W45-W4J. -T45 

T3C 

T3 

WA3-WA31-WB3 

NGDT 

DELTAT 

NG 

NPDT 

DELTAT 

NP 

P41DT 

DELTAT 
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PR TRIM PI... 

LISTING 

5 M68000 

02/19/85 14 5 34 X 16 




s 9? 

P41 

45 

206 

46 0 6 

5666 RQJ 

P45DT 




CNS 

DELTA! 

3 

206 

1794 

S9i 

P45 

56o6 l-:u i 

T3DT 




CN i 

DELTA! 




39 ; 

T3 

*..) o 

206 

4430 

5666 RQJ 

WS3DT 




CN? 

DELTA! 




39? 

WS3 
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PRTRTMPL LISTING 2 M68000 


02/19/85 1 3 5 3 -‘•l X 1 6 


p* 

EXTIME 

CANSTRT. 

MUSTEND 

COMPUTATIONAL 

1 

5666 

0 

5666 

RTTH2-NGG-PCNGG 

-NGDT-NG 


^562 

0 

5666 

P3-WA31-FAR31-H 

-T35-THTA35-W3 

3 

1758 

232 

5666 

PS3-PS3Q2-T3Q2- 

3 

1356 

2016 

5666 

W3 ;|. — P3 1 DT--P3 1 

5 

1273 

3538 

5666 

W XQ 2 - P 3 5 D T P 3 5 

6 

83 0 

0 

3386 

P R3 9 Q 5 “ DH QTH 5 

7 

828 

0 

3120 

PR35Q1-DHQTH3 

8 

598 

1302 

3138 

B1 

? 

598 

230 

3650 

W35C 

10 

390 

2853 

3893 

T0RQ31 

1 1 

366 

568 

3300 

T25Q2-T25-H25 

12 

no C 
OOO 

330 0 

5666 

WS3DT-WS3 

13 

292 

3538 

3138 

WA2 

1 3 

218 

0 

3856 

DELS 

15 

213 

0 

2316 

WF 

16 

160 


5666 

DH35-H39 

J. / 

132 

0 

2530 

H3 

18 

130 

3136 

5330 

WB3 

19 

126 

0 

3300 

H2 


■ACKING SEQUENCE 


' v **** PACKING PROCESSOR 1 

1 , RTTH2-NGC-PCNGC-WA2C-B2-WB25-WA3-TaRQC-NGDT-NG 
MCMCXC3K3K RELATIONAL EQUATIONS ARE 2 

0 ♦ PS3Q2- E: 1 -WA2-H25-H2-H3-T0RGW 1 - WXQ2- WS3DT 
***** PACKING RELATIONAL EQUATIONS 
***** ORDERED EQUATIONS TO BE PACKED ? 

0, PS3Q2“WXQ2-Bl-T0RQ i *l-H25-WS3DT-WA2-H3-H2 
***** NOT PATH ENDPOINT JPS3Q2 
***** ORDERED EQUATIONS TO BE PACKED ♦ 

0, WXQ2-B1-T0RQ31-H25-WS3DT--WA2-H3--H2 
***** TRYING TO INSERT PATH 5 
5* NXQ2-P35DT-P3S 

***** PATH DOES NOT FIT v ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED i 
0* B1-T0RQ31-H25-WS3DT-WA2-H3-H2 
***** TRYING TO INSERT PATH 8 

84 B 1 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED t 
0, TORO^l— H25— WS3DT-NA2— H3--H2 
***** TRYING TO INSERT PATH 10 

10 4 T0RQ31 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED ♦ 

0, H25-WS3DT-WA2-H3-H2 
***** TRYING TO INSERT PATH 11 

11 ♦ H25--T25-T25Q2 

**** PATH DOES NOT FIT » ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED i 
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PRTRTMPL LISTING 


PAGE 


0, W S3 DT —W A2.-H 3 -H2 
***** TRYING TO INSERT PATH 12 

12. WS3DT-WS3 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
rtOK JKMwvC ORDERED EQUATIONS TO BE PACKED X 
0 ♦ WA2-H3-H2 

***** TRYING TO INSERT PATH 13 

13, WA2 

9K3KXC3K3K PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED X 
0 , I-I3--H2 

***** PACKING DELAYED » NOT CRITICAL * H3 
***** ORDERED EQUATIONS TO BE PACKED X 
0. H2. 

***** TRYING TO INSERT PATH 19 
19, 1-12 

***** PATH DOES NOT FIT •/ ADDED TO CARRYOVER* FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED X 
0. H3 

***** TRYING TO INSERT PATH 17 
17, H3 

***** PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 


***** PACKING PROCESSOR 2 

2. P3~WA31-FAR^l-"IHl--T^l--THTA^l-"DH''ll~H' f H--!-K{5"-T‘'15--THTA‘ f 15-N'45™T0RQ' £ l5"-NPDT 


-NP 

***** RELATIONAL 


:quations AREi 


0 , WF“-H3 -DHQTH"I--W X '15C - DH-T5 


3-W‘T 1 -TORCH 1 -HT9- PT 1 DT -P -T5 DT ~W S3 DT ' 


**** CARRYOVER EQUATIONS ARE X 
0 , H3-H2-WA2-WS3DT-H25-TDR0-4 1 -B 1 -WXQ2 
***** PACKING RELATIONAL EQUATIONS 
***** ORDERED EQUATIONS TO BE PACKED X 

0, PS3-W‘ < »l~P^lDT-P‘ e »5DT-DHQTH^-W / l5C-T0RQ^l-WS3DT-WF-DH i J5-H^9-H3 
***** TRYING TO INSERT PATH 3 
3 , PS3-PS3Q2--T3Q2— T3C-T3DT--T3 

2, PS-PSS-WASl—FAR^i— H^l— T^l-THTA i ll-DH‘ ( U.“H / H— H^S—T^S— THTA^S-tkMS—TORQ^S 
--NPDT-NP 

2, F , 3-“PS3" , PS3Q2/"WA31“FAR*ll-H^l“T^l™TITTA^l-"DH^l“H'l^~l--K{5“-T / }5-"THTA^5""N/l5-"TOR0Kt 
--NPDT-NP 

2, P3-“PS3 -"PS3Q2-- WA31“ FAR / 11'-I-K{1~T"11-"THTA^1~DI l'll“' H^PL-HT5”T*15""TI !TA / l5""N' f l5"--T3Q2 
™T0rV(H5-NPDT--NP 

2, P3"PS3"PS3Q2“NA31--FAR‘l:l/~H41~T / Ll™THTA / 41““DH41~-H‘T / i“H ,f t5“T , '}5"-THTA / l5~W / {5"T3Q2 
-T0RCM5-T3C--NPDT--NP 

***** PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED X 

0* W^1HP^1DT-P*5DTH3HQTH4~W45C~T0RQ41-WS3DT-WF-DH45--H i »9--H3 
***** TRYING TO INSERT PATH T 
T , WTl-P^lDT-P^l 

2 , P3—WA3 1 -FAR/1 1 -1-1'! 1 -T'l 1-THT A'T 1 -Dl-K 1 -H - ; H! -H -15 --T H5-T l-l'T M 5- W9 5- N/1 1 -TGRGH5 
-NPDT-NP 

2, P3-WA31-FAR^LI.-l-ITl-THl--TirTATi-DllTl-H4^-IL15-T'15"TTlTA^5~W <: 15--N / ll-T0Ri::K{5 
--NPDT-NP-P41 

***** RATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
'**** ORDERED EQUATIONS TO BE PACKED X 
0 , P' f LI.DT~P / $5DT~-DHQTH / {-~W-T5C“TORG' i 31~WS3DT-WF"-DH'T5—H- < 19~H3 
***** NOT PATH ENDPOINT il :: H IDT 
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PAGE 6 


5!Ok;k>k;k ORDERED EQUATIONS TO BE PACKED i 

()., PT5DT-DHQTHT-WT5C-TQRQ41--WS3DT-WF--DHT5-T-M9-H3 
jkxokjk NOT PATH ENDPOINT JPT5DT 
xxxxx ORDERED EQUATIONS TO BE PACKED i 

()., DHQTH*-WT5C-T0RCm-WS3DT-WF--DHT5“HT9-H3 
TRYING TO INSERT PATH 7 
7 « D H Q T H T - P R T 5 Q 1 

2 , P3-WA3 1 -FART 1 -HI 1 -IT 1-TI-lT AT 1 -DHQTHT--DHT 1 -HII-HI5-TT5-THTAT5-WI5-T0RGI5 
-NPDT-NP 

2* P3-WA31— FART1—HT1—TT1— THTAT1— F :, RT5(3:L~DHQTHT—DHT1~HTT—HT5—TT5— THTAT5-WT5 
-T0RQT5-MPDT-NP 
xxxxx PATH ADDED 

xxxxx UNPACKED PATHS { 3 , T3 ? T , PT 1 > 5 ♦ PIS ? 6 ♦ DHQTH5 ? 8 , B 1 

9 ♦ NTSC ? 1 0 . TORQT 1 i 1 1 ♦ H25 5 12 * WS3 LI. 3, WA2 ? IT ♦ DEL 2 % 15 ♦ WF * 16 , 1-IT9 
17 ♦ H3 5 18 * WB3 * 19 ♦ H2 

xxxxx ORDERED EQUATIONS TO BE PACKED t 
0« WI5C~TOROW1~WS3DT-WF“DHT5-HI9-H3 
xxxxx TRYING TO INSERT PATH 9 
9* WT5C 

xxxxx PATH DOES NOT FIT? ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED l 
0 * TORQT 1 -WS3DT~WF-DHT5~HI9~H3 
xxxxx TRYING TO INSERT PATH 10 
10* TORQT 1 

xxxxx PATH DOES NOT FIT? ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED t 
(! * WS3DT-WF-DHI5-HT9-II3 
xxxxx TRYING TO INSERT PATH 12 
12* WS3DT--WS3 

2, P3-- W A3 1 -FART 1--HT 1 -TT 1 -THT AT 1 -PRTSQ 1 -DHQTHT-DHT 1 --HII— HIS— TI5-THT AI5- WI5 
— T0RQI5-NPDT— WS3DT-NP 

xxxxx PATH DOES NOT FIT? ADDED TO CARRYOVER F"OR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED t 
0 * WF-DHT5— HT9-H3 
xxxxx TRYING TO INSERT PATH 15 
15* WF 

2 , P3-WA3 1 -WF-FART 1 -HI 1 -TT 1 --THTAT 1 -PRTSQl-DiTQTHT- DHT 1-HTT-HT5-TT5- T H f ATS 
-WT5-T0RQT5-NPDT-NP 
xxxxx PATH ADDED 

xxxxx UNPACKED PATHS t 3 ♦ T3 ? T * PT1 ? 5 <■ PT5 ? 6 * DITQTH5 5 8 * B1 

9 4 WT5C ? 1 0 * TORQT 1 ? 1 1 . H2S ? 1 2 . WS3 f 13 « WA2 Li. T * DEI- 2 ? 1 6 * IT T9 f 1 7 * IT3 
18.WB3L19 *H2 

xxxxx ORDERED EQUATIONS TO BE PACKED i 
0 ♦ DHT5-HT9-H3 

xxxxx TRYING TO INSERT PATH 16 
16* DHI5-HI9 

2* P3-WA81-WF-FART1-HT1-TT1-THTAT:I,-PRT5Q1--DHQTHT-DHT1-HTT-HT5-TT5-THTAT5 
-WT5-T0RQT5-NPDT-NP-HT9 

xxxxx PATH DOES NOT FIT? ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED 5 
0 * HT9-IT3 

xxxxx TRYING TO INSERT PATH 16 
16, HI9-DHI5 

2 , P3-WA3 1 -WF -FART 1 -HT 1-TI 1 -THT AT 1 -PRT5Q 1 -DHQTHT-DHT 1 -HTT-HT5-TT5-THTATS 
—WI5— TGRQI5— NPDT— NP— HI9 

xxxxx PATH DOES NOT FIT? ADDED TO CARRYOVER FOR NEXT PROCESSOR 
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PAGE 


xokjkkcjk ORDERED EQUATIONS TO BE PACKED ? 

0, H3 

kxxx TRYING TO INSERT PATH 17 
17. H3 

xxxxx PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** PACKING CARRYOVER EQUATIONS 
***** ORDERED EQUATIONS TO BE PACKED t 
0. WXQ2--B 1 --TORCH 1 -H25--NS3DT-NA2-I I3--H2 
xxxxx TRYING TO INSERT PATH 5 
5 . P'15--PT5DT--HX.Q2 

xxxxx PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED t 
0 4 B1-T0RQ < H~H25-WS3DT-WA2-H3~H2 
9KMC3K3K9K TRYING TO INSERT PATH 8 
8. 81 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED S 
0 4 TORCH 1 -H25-WS3DT-WA2-H3-H2 
***** TRYING TO INSERT PATH 10 
10. TORCH 1 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED t 
0. H25-WS3DT-WA2-H3-H2 
xxxxx TRYING TO INSERT PATH 11 
1 1 4 H25-T25-T25Q2 

xxxxx PATH . DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
iiox:toi;;K ORDERED EQUATIONS TO BE PACKED t 
0 4 WS3DT-WA2-H3-H2 
xxxx TRYING TO INSERT PATH 12 
12 * MS3-MS3DT 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED t 
0. WA2-H3 -H2 

xxxxx TRYING TO INSERT PATH 13 
13. WA2 

xxxxx PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED $ 

0. H3--H2 

xxxxx TRYING TO INSERT PATH 17 
I/. H3 

xxxxx PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR 
xxxxx ORDERED EQUATIONS TO BE PACKED;: 

0, H2 

xxxxx TRYING TO INSERT PATH 19 
19* IT 2 

:5ox!K>xik PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 


xxxxx PACKING PROCESSOR 3 

3. FS3-PS3Q2-T3Q2-T3C-T3DT-T3 
RELATIONAL EQUATIONS ARE t 
0* T25Q2 

V/ V/ VJJ V/ CARRYOVER EQUATIONS ARE i 

0 4 H2-H3-WA2-WS3DT--H25-T0RQ^l-Bl-WXQ2-H3-H^9--DH / l5-WS3DT-TaRQ' < H-W^ 
-PS 3 

PACKING REI._ATT.ONAL EQUATIONS 
xxxxx ORDERED EQUATIONS TO BE PACKED i 


SC-PHl 


19 
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0, T25Q2 

xxxxx TRYING TO INSERT PATH 11 
11 , T25Q2-T25-H25 

3» PS3-PS3Q2-T25Q2-T3Q2-T3C-T3DT-T3 
3 * PS3— PS3Q2—T25Q2—T 25—T3G2 -T3C-T3DT ~T 3 
3 4 PS3-PS3Q2-T25Q2-T25-H25-T3Q2~T3C-T3DT-T3 
***** PATH ADDED 

***** l JNPACKED PATHS X *1 •> P4 1 5 5 •» P45 ? 6 , DHQTH5 ? 8 ♦ B 1 * 9 ♦ WTbC 

1 o » TORCH 1 5 12 * WS3 J 13 ♦ WA2 5 14 ♦ DEL 2 5 16 > H49 1 17 , H3 $ 18 « WB3 5 19 .H2 
***** PACKING CARRYOVER EQUATIONS 
xxxx* ORDERED EQUATIONS TO BE PACKED t 

0 4 WT1-WXQ2H31H4T5C-TORQT1-WS3DTH4A2H3H45-H49--H3--H2 
xxxxiK TRYING TO INSERT PATH T 

1, PTl-PTlDT-WTl 

3* PS3~PS3Q2-T25Q2-T25~H25~T3Q2-T3C-T3DT-T3-P i U 

3 4 PS3--PS3Q2 -T25Q2— T25“"H25--T3Q2“-T3C"T3DT--T3“P / 1 1 DT-P4 1 

3 , PS3 ••••PS3Q2-T25Q2-T25 -H25"-T3Q2""T3C”T3DT--T3 -•W / LL™P-^1DT“P / !1 
xxxx* PATH ADDED 
xxxxx RELATIONAL EQUATIONS ARE t 

0 4 TDRQ41-P45DT 

3K3KX3KMC UNPACKED PATHS: 5.P45;64DHQTH3?8.B1 19* W45C; 10 «T0RGM1 
1 2 » WS3 f 1 3 4 WA2 5 1 4 * DEL2 516* I I T9 P 17, H3 f 1 8 * UB3 > 1 9 4 H2 
xxxxx ORDERED EQUATIONS TO BE PACKED: 

0 * PT5DT--T0RQT1 
***** NOT PATH ENDPOINT JP45DT 
***** ORDERED EQUATIONS , TO BE PACKED : 

0 4 TORCH! 

xxxxx TRYING TO INSERT PATH 10 
10 4 T0RQ41 ' ’ 

3 .> PS3-PS3Q2-T25Q2-T25HH23-T3a2-T3C-T3DT-T3-W41-T0RQ41~P41DT--P41 
***** PATH ADDED 

***** UNPACKED PATHS : 5 ♦ P45 > 6 •> DHQTH5 5 8 4 B1 5 9 4 WT5C 5 12 ♦ MS3 
13 4 WAS t IT 4 DEL.2 ? 16 ♦ I-IT9 5 17 ♦ H3 5 18 ♦ WB3 ? 19 4 H2 
***** ORDERED EQUATIONS TO BE PACKED: 

(), WXQ2-B1-W45C-WS3DT-WA2-DH45-H49-H3-H2 
***** TRYING TO INSERT PATH 5 
54 P45-P45DT-WXQ2 

PS3-PS3Q2-T25Q2-T25-H25-T3Q2-T3C-T3DT-T3-WT1-T0RQT1-PT1DT--PT1-PT5 
PS3™PS3R2~T25R2— T25—H25~T3Q2“T3C— T3DT“T3-“WT1~T0RQT1— PT1DT— PT1“PT5DT 

--PT5 

PS3-PS3Q2-T25Q2'-T25-H25“T3Q2“T3C-T3DT-T3-WT1-T0RQT1-PT1DT-PT1-WXQ2“PT5DT 
-PT5 

***** PATH ADDED 

xxxxx RELATIONAL EQUATIONS ARE : 

0 * WA2-WB3 

xxxxx UNPACKED PATHS : 6 •> DHQTH5 1 8 * B 1 % 9 ♦ WT5C 4 124 WS3 5 13 ♦ WA2 
1 T 4 DEL 2 ; 16 4 HT9 ? 1 7. H3 ? 18 4 WB3 519 4 H2 
xxxxx ORDERED EQUATIONS TO BE PACKED : 

0 4 WA2-WB3 

xxxxx PACKING DELAYED ■> NOT CRITICAL: WA2 
***** ORDERED EQUATIONS TO BE PACKED: 

0* WB3 

***** TRYING TO INSERT PATH 18 
18 4 WB3 

3 4 PS3-PS3Q2— T25Q2— T25— H25~T3Q2~T3C— T3DT—T3— WT1— TORQT1— F’TIDT— PT1—WXQ2—WB3 
—PT5DT—P45 


O 4 

o A 


3. 
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***** PATH ADDED 

***** RELATIONAL EQUATIONS ARE ? 

0* WA2-WS3DT 

***** UNPACKED PATHS ? 6 , DHQTH5 ? 8 * B1 * 9 * W45C f 12 ♦ WS3 513* WA2 
1 T ♦ DEL2 * 1 6 ♦ FH-9 1 17* H3 * 1 9 * H2 
***** ORDERED EQUATIONS TO BE PACKED? 

()« WS3DT-WA2 

***** TRYING TO INSERT PATH 12 
12 4 WS3DT-WS3 

3* pS3“-PS3Q2“T25Q2"T25--H25-"T3Q2“T3C--T3DT""T3“W / J:!.-"T0RCV}l~P'HlDT“P'tl"-WXQ2--WB3 
-PT5DT--W33DT--P-T5 

***** PATH DOES NOT FIT v ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED? 

0 * WA2 

***** PACKING DELAYED * NOT CRITICAL? WA2 
***** ORDERED EQUATIONS TO BE PACKED ? 

0* Bl-W45C-WS3DT-WA2-DH / l5-“H i 19“-H3-H2 
***** TRYING TO INSERT PATH 8 
8* 131 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED? 

0 , WT5C“ NS3DT“WA2“DH / f5-'HT9“H3-H2 
***** TRYING TO INSERT PATH 9 
9, UH5C 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED? 

0 ♦ WS3DT—WA2— DH^S—ITW— H3--H2 
' v **** TRYING TO INSERT PATH 12 

1 2 ■> MS3-MS3DT 

3, PS3-PS3Q2-T25Q2-T25-H25-T3Q2-T3C-T3DT~T3-W / H-T0RQ^1-P41DT-P41-WXQ2-WB3 
""PTSDT-WSSDT-P^S 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR ‘ 

***** ORDERED EQUATIONS TO BE PACKED? 

0 4 HA2-DHH5-I-R9 - H3-H2 
XX*X* TRYING TO INSERT PATH 13 

1 3 •> WA2 

***** PATH DOES NOT FIT* ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED ? 

0 * DI-H5--HH9- H3--H2 
***** TRYING TO INSERT PATH 16 
16 4 HT9--DHH5 

PS3H*S3Q2-T25Q2-T25-H25-T3Q2-T3C-T3DT-T3“W / »l-T0RQ41~P41DT-P‘U-WXa2~WB3 

-P^SDT-PHS-HHS 

***** PATH DOES NOT FIT * ADDED TO CARRYOVER FOR NEXT PROCESSOR. 

***** ORDERED EQUATIONS TO BE PACKED? 

0 4 HT9--H3-H2 

***** TRYING TO INSERT PATH 16 


16, l-K 


DI-IT 5 


3, PS3"-PS3Q2-T25Q2“T25-H25~T3a2-T3C-T3DT-T3-W / H~T0RQ < H~P‘‘llDT-P^l-WXQ2-WB3 
- F’ T 5 D T - F’ '•} 5 - H T 9 

***** PATH DOES NOT FIT » ADDED TO CARRYOVER FOR NEXT PROCESSOR 
***** ORDERED EQUATIONS TO BE PACKED? 

0 > H3--H2 
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~ W B 3— P A 5 D T — P A 5 
ioox PATH ADDED 

**** UNPACKED PATHS t 6 , DHQTH5 5 8 * B 1 5 9 , WA5C 5 12 ♦ NS 3 513 . WA2 
1 A * DEL2 5 16 » I-IA9 5 19 . 1-12 
***** ORDERED EQUATIONS TO BE PACKED t 


0* H2 

***** TRYING TO INSERT PATH 19 
19 4 H2 

3., H3H-I2 ••-PS3~PS3Q2—T25Q2~T25--H25— T3Q2— T3C—T3DT-T3"-WA1— TGRQA1-PA1DT--PA1 
-NX Q2-WB3-P AS DT -P AS 
***** PATH ADDED 

***** UNPACKED PATHS t 6 * DHQTH5 5 3 ♦ B1 5 9 , WA5C 5 12 , WS3 5 13 * HA2 
1A , DEL 25 16 ♦HAS 1 


***** PACKING PROCESSOR A 
A < PRA9Q5-DHQTH5 
***** RELATIONAL EQUATIONS ARE i 
0 , WA5C-DHAS 

***** CARRYC )VER EQUATIONS ARE: 

0, HA9--DHA5--WA2-WS3DT-WA5C-B1-WS3DT 
***** PACKING RELATIONAL EQUATIONS 
***** ORDERED EQUATIONS TO BE PACKED t 
0 * WA5C-DHA5 

***** TRYING TO INSERT PATH 9 


94 HA5C 

A 4 PRA9Q5-WA5C-DHQTH5 


***** PATH ADDED 

***** UNPACKED PATHS X 8 ♦ B1 5 12 ♦ WS3 5 13 ♦ MA2 X 1A 4 DEL2 5 16 4 1-IA9 
JKMC3KX ORDERED EQUATIONS TO-BE PACKED: 

0 4 DI-IAS 

***** TRYING TO INSERT PATH 16 


16, DHA5-HA9 

A 4 PRA9Q5-NA5C--DHQTH5--DHA5 
A 4 PRA9Q5-WA5C-DHQTH5-DHA5-HA9 
***** PATH ADDED 

***** UNPACKED PATHS : 8 4 B1 5 12 4 WS3 5 13 ♦ WA2 5 1 A « DEL2 
***** PACKING CARRYOVER EQUATIONS 
***** ORDERED EQUATIONS TO BE PACKED: 

0, B 1 -WS3DT-WA2 
***** TRYING TO INSERT PATH 8 
8 ♦ B1 

A 4 PRA9Q5--.WA5C-DHQTH5--B1 --DHA5--HA9 
***** PATH ADDED 

***** UNPACKED PATHS : 12 , WS3 5 13, WA2 51 A * DEI..2 
***** ORDERED EQUATIONS TO BE PACKED: 

0 , WS3DT-WA2 

***** TRYING TO INSERT PATH 12 


12, WS3-WS3DT 

A * PRA9Q5--WA5C ~DHQTH5--B1-DHA5--HA9--W33 
A, PRA9Q5—WA5C— DHQTH5—B1—DHA5— HA9—WS3DT—WS3 
***** PATH ADDED 

***** UNPACKED PATHS : 1 3 , WA2 5 1 A 4 DEL2 
***** ORDERED EQUATIONS TO BE PACKED : 

0, WA2. 

.•!■£ )!\ y<-. /!{ ' TRYING TO INSERT PATH 13 

13, WA2 
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PRT9Q5-WT5GH3HQTH5-B1-DHT5-WA2-H49-WS3DT-WS3 
xxxxx PATH ADDED 
xxxx RELATIONAL. EQUATIONS ARE t 
0, DEL 2 

xxxxx UNPACKED PATHS t 1T.DEL2 
xxxxx ORDERED EQUATIONS TO BE PACKED i 
0, DEL 2 

xxxxx TRYING TO INSERT PATH IT 
IT, DEL 2 

T , PRT9Q5-WT5C-DHQTH5-B1-DEL2-DHT5-WA2-HT9--WS3DT-WS3 
xxxxx PATH ADDED 
xxxxx UNPACKED PATHS t NONE 

:t< SMALLEST PATH EXECUTION TIMES 3 (HZ? UNUSED TIME ON OTHER PATHS* 186 
PACKING SEQUENCE COMPLETE 
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RESULT 

PROCESS 

CALC 

OR EQUAT 
START 

ION — 
END 

MAX 

RTTH2 

770 

0 

770 

.... . TT -, tf .- 4 . 

770 

NGC 



434 

770 

1204 

1204 

PCNGC 

198 

........ — i., 

1204 

1402 

1402 

WA2C 

— • i 

2136 

1402 

3538 
**** **** **** **“ 

3538 

82 

610 

3338 

4148 

4148 

WB25 

... . . „ » 

122 

i. , 

4148 

4270 

4270 

WAS 

“ T 

30 

„ HM A. . 

4270 

4300 

4300 

TORQC 

594 

*130 0 

... .... .» 

4894 
._ .... . 

4894 



r 

NGDT 

..... ,1. . 

566 

r"" 

4894 

T 

5460 

5460 

NG 



— “"•—V 

2.06 

5460 

5666 

5666 


< SENT TO > < NEEDS - > 

PR EQUATION REQ PR EQUATION AVAIL 


4 WA2 

— +- 

3706 


— 



4 81 


1438 




3 WXQ2 


4138 

3 PS3Q2 


594 

4 WA2 ' 

— i — 

3706 












4 B1 


2036 




4 WA2 


3998 

4 WS3DT 


48 66 

4 WA2 


3998 




3 1-13 
3 H2 

142 

268 





3 1125 


960 




4 WA2 


3998 


-4 — 


3 TOR 0.41 


3620 


x 

x 

X 


X SC SC SC SC SC IK IK IK SC 

PROCESSOR 2 
UTILIZATION : 98% 

SC X X SC SC SC SC X SC X 


X 

sc 

X 

sc 


RESULT 

- PROCESS 
CALC 

OR EQUAL 
START 

ION 

END 

MAX 

P3 

242 

0 

2.42 

304 

WA31 

970 

242 

1212 

1274 

WF 

214 

1212 

1426 

1488 

FAR41 

224 

1426 

1650 

1712 

1-141 

336 

1650 

1986 

2048 

v ‘4 1 

122 

1986 

2108 

2170 

THTA41 


2108 

2230 

2292 


< SENT TO > < NEEDS -> 

PR EQUATION REQ PR EQUATION AVAIL 


3 PS 3 


268 


3 P41DT 

4 WS3DT 


3620 

4866 













3 H3 142 

3 P410T 


3620 


3 W41 


2392 
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. LIS 

TING : M630D0 

02/1 
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PAGE 
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PRT5Q1 

230 

2230 2 TOO 

2522 





DHQTHT 

598 

2T60 3050 

3120 






di-iti 

118 

3058 3176 

3238 

3 TORGT1 

3230 


i — 


i-itt 

5-1 

3176 3230 

3292 





HT5 

11T 

3230 33TT 

3T06 

T HT9 

3998 




TT5 

122 

33 TT 3T66 

3523 

3 PT5DT 

T866 




TIITAT5 

122 

3066 3588 

3650 

T- DIMS 

3588 




NTS 

85 T 

3583 TTT2 

T50T 

3 PT5DT 

T866 

T WT5C 


823 

TGRQT5 

390 

TOT 2 0832 

T89T 



T DI-IT 5 


3706 

NPDT 

566 

T832 5398 

5T60 






NP 

206 

5393 56 OT 

5666 







:i< sc :x >t< :« s< sc xc :k :k ;k s< 
i« PROCESSOR 3 * 

is UTILIZATION. J 97% x< 
;k ;k :k ;k ;k ;k x ;k »: :k s< k 


. PROCESSOR EQUATION — > 


RESULT 

CALC 

START 

END 

MAX 

H3 

1T2 

0 

1T2 

266 

H2 

126 

1T2 

268 

392 

PS3 

122 

268 

390 

51 T 

PS3Q2 - 

20 T 

390 

59 T 

713 

T25Q2 

13T 

59T 

728 

352 

1 X..J 

118 

723 

8T6 

970 

1-125 

1 IT 

8T6 

960 

108T 

T3Q2 

598 

960 

1558 

1682 

T3C 

13T 

1553 

1692 

1816 

T3DT 

T9T 

1692 

2186 

2310 

To 

206 

2186 

2392 

2516 

/•} 1 

838 

2392 

3230 

335T 




< SENT TO > < NEEDS > 

PR EQUATION REQ PR EQUATION AVAIL 


1 TORQC 

2 l-ITl 


T300 

1650 




1 TORQC 


T300 

2 P3 






2T2 

1 WA2C 


IT 02 






*“* f.— ■•“**“*{ **"**“.“-*““**”•— — .f— ■ 


1 TORQC T300 





i 

t 

t 

j 

i 

! 

i 

+ 

i 

i 

1 

! 

2 
i 
i 
i 












—■■i 

2 TI-ITATl 

2230 
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TING X H 
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PAGE 

10 


TORSO 1 390 

3230 

3620 

3700 

1 NGDT 

0890 

2 DM0 1 


3176 

01 DT 3:12 

3620 

3932 

0056 



2 WA31 


1212 







2 TOl 


2103 

POl 206 

3932 

0133 

0262 






WXQ.2 598 

0133 

0786 

0736 

0860 



1 WA2C 


3533 

WB3 130 

0866 

0990 

0 WS3DT 

0866 

0 WA2 


3998 

P05DT '17 0 

0866 

5336 

5060 



2 T05 


3066 







2 W05 


/J./. j /}.•”> 







0 WA2 


3993 

P'05 206 

5336 

5502 

5666 






;k ;;< :k ik n< ik jk ik 

ik IK IK IK 








* PROCESSOR 0 

IK 








;k UTILIZATION J 53% sc 








»: x< ik ix ik x< :k ik 

IK IK IK XC 








DDnrnroci 

n rrm iatthm 

... 

< SENT TO 



i:r r\o 


*•. I KUULOO 

t-.uewn i . 

Jk.UIX 


/ 

*■. nil".! 

L 1 ./ 


RESULT CALC 

START 

END 

MAX 

PR EQUATION 

REQ 

PR EQUATION 

AVAIL 

PR09G5 230 

0 

230 

1710 




— 


, i'05C 598 

230 

828 

2312 

2 W05 

3538 



DHQTH5 610 

823 

1038 

2922 






01 ’ 598 

1038 

2036 

3520 

1 WB25 

0108 

1 PCNGC 


1002 

DEL2 218 

2036 

2250 

3738 






DILI 5 118 

3583 

3706 

3356 

2 TORSOS 

0002 

2 THTA05 


3583 

WA2 292 

3706 

3998 

0108 

1 WB25 

0108 

1 RTTH2 


770 





1 WAS 

0270 

1 WA2C 







1 TORQC 

0300 








3 WB3 

0736 








3 P05DT 

0866 




HO 9 02 

3998 

oooo 

5330 



2 H05 • 


3300 

WS3DT 130 

0366 

0996 

5060 



1 WAS 


0300 







2 WA31 


1212 







3 WB3 


0866 

WS3 206 

0996 

5202 

5666 
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TABLE I. - RTMPL - GENERATED INFORMATION 


Label 

Operation 

Result 

Arguments 


S$12 



z 

8 



■ 3 |: 

Ri,y 

16 


STORE 

mm 

R1 

8 


TABLE II. - PARTITIONING INFORMATION 


Equation 

label 

Dependent 

arguments 

Independent 

arguments 

1 

m 

1 

■ 

Must 

end 3 

X ! 

(none) 

u. x 5 

32 

0 

32 

64 

X 2 

(none) 

X 5 

64 

0 

64 

64 

X 3 

X 2 

(none) 

32 

64 

96 

112 

X 4 

CNJ 

X 

(none) 

48 

64 

112 

112 

X 5 

* 3 . X 4 

X 5 

48 

112 

160 

160 


information used for packing only. 


TABLE III. - PACKING ALGORITHM RESULTS FOR TURBOJET MODEL 


B 

Processors 

required 

Processor percent utilization 

P number 1 

P number 2 

P number 3 

P number 4 

5666 

4 

100 

98 

97 

53 

10000 

2 

99 

98 

— 

— 

19568 

1 

100 

— 

— 

— 
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1 Dual generator simulation 1 

1 Computer #1 Computer #2 1 

it 

Sample u 

Sample u 


l A , T 

Compute G1 

Compute G2 

j Controller 

-L- 

■ 

Output G1 

Output G 2 

|Gi. G 2 


(a) Decoupled computational paths. 


Computer h 


Aircraft simulation 

Computer #2 


5AT Read E(iAT ) \ 
Compute A ((i+5)AT) 


Read A (iAT), . 

Compute E ((i+l)AT) 

Use A (iAT)/ v 

Compute E((i+2)AT) 


I 


Read E ((i+5) At) 
Compute A ((i+10)AT; 


Use A (d+5)AT) 
Compute E ((i+7)ATj 



(b) Loosely coupled computational paths. 
Figure 1.- Physical parallelism. 





Computer ^2 


Computer #1 


Computer #3 


X U 

Compute X 2 

X 2 , 


. Sample u 
Compute Xj 


Wait Xi 
Compute X 4 
Wait X 3 

*3 



Wait X 2 
Compute X 3 



Compute X 5 




(a) Closely coupled paths. 


Computer ^1 
Computer ^2 
Computer #3 

Computer h 
Computer ^2 

Figure 2. - Partitioning and packing closely coupled equations. 


Time — 

TX2 

TX4 

TX5 

T 1 

i 

(Idle) 

(Wait) 

T X3 1 

ill (Idle) 


f (Packing) 



(b) Packing. 









Figure 3. - Path identification Algorithm. 














*See explanation in test. 


Figure 4. - Packing Algorithm. 














Figure 4.- Continued. 
















■"See explanation in text 


Figure 5. - Test fit Algorithm. 
















CANSTART (a) 


CANEND (a) 
MUSTEND (b) 


CANSTART (b) 
CANEND (b) 


MUSTEND (b) 



CANSTART (0 


CANEND (0 
MUSTEND (0 


CANSTART (a) 


CANEND (a) 
MUSTEND (a) 
CANSTART (0 


J” 


CANEND (0 
MUSTEND (0 
CANSTART (b) 


CANEND (b) 
MUSTEND (b) 


Processor before inserting 
equation 


Equation to be inserted 


Processor after insertion 


Figure 6. - The affect of insertion on equation attributes 
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